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Entanglement renormalization is a real-space renormalization group (RG) transformation for 
quantum many-body systems. It generates the multi-scale entanglement renormalization ansatz 
(MERA), a tensor network capable of efficiently describing a large class of many-body ground 
states, including those of systems at a quantum critical point or with topological order. The MERA 
has also been proposed to be a discrete realization of the holographic principle of string theory. In 
this paper we propose the use of symmetric tensors as a mechanism to build a symmetry protected 
RG flow, and discuss two important applications of this construction. First, we argue that symmetry 
protected entanglement renormalization produces the proper structure of RG fixed-points, namely a 
fixed-point for each symmetry protected phase. Second, in the context of holography, we show that 
by using symmetric tensors, a global symmetry at the boundary becomes a local symmetry in the 
bulk, thus explicitly realizing in the MERA a characteristic feature of the AdS/CFT correspondence. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Hk 



Renormalization [l[ is fundamental to our understand- 
ing of quantum many-body systems. The renormaliza- 
tion group (RG) explores how the behavior of an ex- 
tended system depends on the scale of observation. Im- 
portant concepts such as universality, criticality or stabil- 
ity of phases are then explained in terms of the existence 
of fixed-points of the RG flow. On the other hand, in 
holographic string theory constructions 0, Q , an addi- 
tional dimension corresponding to scale (or RG flow) is 
used to regard a many-body system as the boundary dual 
of a higher dimensional, bulk theory. 

Both renormalization and holography can be realized 
non-pcrturbativcly, on the lattice, using entanglement 
renormalization, a real-space coarse-graining transforma- 
tion for many-body wave- functions 0, S|. A key as- 
pect of entanglement renormalization is the removal of 
short-range entanglement using so-called disentanglers 
Fig.QJa). This leads to an efficient description of many- 
body states in terms of the multi-scale entanglement 
renormalization ansatz (MERA) f3], a tensor network 
that spans an additional dimension corresponding to 
scale, Fig. [T](b). The MERA has been applied to the ex- 
ploration of frustrated antiferromagnets 0], interacting 
fermio ns B , topological order [j| [lij] , quantum critical- 
ity (Til flij and, more recently, holography [l3|, 14 [. 

In this paper we propose the use of symmetric tensors 
as a way to protect a global symmetry during the RG 
flow, and present two applications of this construction in 
the context of renormalization and holography. Specifi- 
cally, we will explain how to use symmetric tensors to: 
(i) obtain the proper structure of RG fixed-points corre- 
sponding to symmetry-protected phases; and (ii) obtain 
a holographic description where a global symmetry at the 
boundary turns into a local symmetry in the bulk. 




FIG. I: (a) Coarse-graining transformation, made of dis- 
entanglers u and isometries w, that maps a lattice L into 
a coarse-grained lattice £ ' . (b) The MERA is obtained by 
repeatedly applying the coarse-graining transformation, (c) 
Constraints satisfied by CJ-symmetric disentanglers and isome- 
tries, Eq. 01 A global symmetry Q can be enforced by using 
C/-symmetric tensors: equality (b) = (d) follows from the sym- 
metry constraints in (c); (d) = (e) is obtained using V„ V g = I 
on each bond index, which only leaves V g acting on the open 
indices; therefore (b) = (e), which amounts to Eq. [2] 



Symmetry protected entanglement renormalization. — 
Let £ be a lattice made of L sites, each described by a 
vector space V, and let V g : V — > V be a (unitary or 
anti- unitary) linear representation of the group Q, 

V f V g = V h ^f-g = h, V/,.ge0. (1) 

We say that a many-body wave- function \^f) G V® L on 
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the lattice C is invariant under the global symmetry Q, 
or Q -symmetric (IBj . if 



(U g )® L \#) = \*) 



(2) 



Similarly, we say that a tensor is ^-symmetric if it has 
(unitary or anti-unitary) linear representations of Q on 
all of its indices, and is invariant under its action. In the 
MERA [6], where the disentanglers u and isometries w, 
see Fig. QJa), implement maps between one or two sites, 
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invariance under Q is expressed as (see also Fig. [ljc)) 
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The result of coarse-graining a ^/-symmetric wave- 
function |W) using ^-symmetric disentanglers u sym and 
isometries W sym is also a ^-symmetric state, and we thus 
obtain a symmetry protected RG flow. In particular, 
Fig.QJb,d,e) shows that the MERA, which is built by con- 
catenating several coarse-graining transformations, rep- 
resents then a (/-symmetric state |\&) G V® L [lq . 

Application I: symmetry protected RG flow. — For con- 
creteness, let us specialize to time reversal symmetry, 
Q = Z% , acting on a lattice £ where each site is described 
by two spin-1/2 degrees of freedom, V = (C2)/ <8> (C2) r - 
On each spin, time reversal acts anti-unitarily by means 
of operator T, 



T - 



i&ylC, 



(6) 



where a y is a Pauli matrix and AC denotes complex con- 
jugation. Since T"T — —I, time reversal acts on a sin- 
gle spin-i as a projective representation of Z2. How- 
ever, when acting simultaneously on two spins, we have 
(f ® 2 )(f ® 2 ) = (-i)® 2 = (i)® 2 , and therefore time rever- 
sal acts on each site of lattice C by means of a linear rep- 
resentation of Z2. A singlet state \^p~) = ^(|01) — 1 10)) 
of two spins is invariant under time reversal, as follows 
from (ia y )® 2 \ip~) = \ip~). We then consider the follow- 
ing two Z^-symmetric many-body states, Fig.[2Ia,b), 



|* prod ) EE C ./-,. 

sec 

I^aklt) = \lp~)s,r;a+l,V 

sec 



(7) 
(8) 



In the product state l^prod) each site s of C is in a sin- 
glet state and therefore there is no entanglement 
between different sites. Instead, in the entangled state 
$ AKLT ), which is related to the AKLT state 17j , a singlet 
state entangles the right spin of site s with the left spin 



^ ^ f"'*"^ — 1 4 1 ) ^ ^ ^prod 

) = i 



(c) 



(e) ^ (f) 





FIG. 2: (a) Pictorical representation of the singlet state 
4>~ (in bra and ket forms) and of the identity operator I on 
the space C2 of a spin-i degree of freedom, (b) States ^ pPO d 
and $ aklt j Eqs. 17181 made of singlet states ip~ . (c) Zj- 
symmetric disentanglers and isometries. Notice that fi Bym is 
the identity operator (it cannot remove short-range entangle- 
ment, which is protected by the symmetry), (d) As a result, 
the coarse-graining transformation only succeeds at removing 
every second singlet state In red are examples of a sin- 

glet state that cannot be removed (left) and one that can be 
removed (right) during coarse- graining, (e) Non-symmetric 
disentanglers and isometries. By allowing some of their in- 
dices to contain projective representations of Zf* , instead of 
linear representations, disentanglers can now eliminate short- 
range entanglement, (f) As a result, the state 4/aklt becomes 
the state $ prod after just one coarse-graining transformation, 
(g) The ground state *&aklt is thus a fixed-point of the RG 
flow with symmetry protection. 



of site s + 1 , and therefore there is entanglement between 
nearest neighbor sites. It is easy to see that both wave- 
functions are unique ground states of Z^-symmetric, lo- 
cal, gapped Hamiltonians, namely H plod = J^sec &s,l'&s,r 
and H ak.lt = 12 sec &s,r'<?s+i,i- [Indeed, each interaction 
a ■ a acts on a different pair of spins, and \ip~) is the 
unique ground state of a ■ a] . 
States * prod ) and * 

aklt) arc interesting because in a 
system invariant under time inversion they belong to two 
different phases 18, ljjj . with I^aklt) corresponding to 
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the Haldane phase 17 1. Following Refs. Tj^- 2(|, we say 
that two local, gapped Hamiltonians correspond to the 
same phase of matter if we can smoothly deform one into 
the other by a means of a path of local, gapped Hamil- 
tonians. Similarly, we say that two ground states are in 
the same phase if they are the unique ground states of 
those Hamiltonians. In absence of symmetry, in one spa- 
tial dimension there is only one phase }18l420| | , and we can 
choose the product ground state l^p^d) as its representa- 
tive. However, when we add the restriction that a global, 
on-site symmetry Q must be preserved at all times, then 
one obtains several inequivalent phases, called symme- 
try protected topological (SPT) phases, classified by the 
second group cohomology H 2 {Q,\](1)) of the symmetry 
group G 18-2(|- For time reversal symmetry there are 
only two SPT phases, because the second group coho- 
mology of Zj is Z 2 , i? 2 (Z^,U(l)) = Z 2 , and |* prod ) and 
I^aklt) are their representatives. 

Our goal is to use symmetry protected entanglement 
renormalization to reproduce the classification of SPT 
phases from the perspective of the RG flow. We have 
found that, indeed, states |\l/ prod } and I^aklt) correspond 
to two inequivalent fixed points of the RG flow when 
Z^-symmetric disentanglers and isometries are used in 
the coarse-graining transformation; and that they both 
flow to the same fixed-point if non-symmetric tensors 
are used. The RG flow of l^p^d) is trivial: a product 
state is transformed into a product state [2^, and we 
will not discuss it here. Let us focus instead on state 
Y&aklt)- Fig. [2Uc),(d) shows a specific choice of Zj- 
symmetric disentangler and isometry (u aym ,u> aym ) that 
coarse-grain state \^aklt) into a state locally identical 
to \^aklt)- In other words, \^aklt) is a fixed-point 
of the symmetry protected RG flow, Fig. EJg). In con- 
trast, when non-symmetric disentanglers and isometries 
are allowed, we can find a pair (u, w) that coarse-grains 
Y&aklt) into the product state IvtVod) after just one step 
of coarse-graining, Fig. |2je),(f). 

It is important to emphasize that although the ten- 
sors (it, w) in Fig. [5Je) are not Z^-symmetric (some of 
their indices carry projective representations of Z^), the 
resulting coarse-graining transformation in Fig. [2](f ) is 
fully Z^-symmetric. In other words, it is only when we 
attempt to implement the Z^-symmetric coarse-graining 
transformation in Fig.[2]Jf) by means of local transforma- 
tion (disentanglers and isometries) that we are forced to 
use projective representations and thus non-symmetric 
tensors. 

A byproduct of the above analysis is an exact MERA 
representation of I^aklt) in terms of the pair of fixed- 
point, Z^-symmetric tensors (w aym ,w aym ). This construc- 
tion can be extended to an arbitrary group Q, by replac- 
ing the singlet state of two spin-i degrees of free- 
dom with an entangled state of two (generalized) spins 
that is invariant under the action of the suitable projec- 
tive representations of Q acting on the two spins, and 




FIG. 3: (a) MERA for the state |*> of a lattice L with 
L = 32 sites in d = 1 space dimensions. This holographic 
tensor network expands one additional dimension, with the 
radial direction corresponding to scale. The MERA can be 
used to represent the ground state of a quantum critical sys- 
tem, corresponding to a CFT, and its network of tensors is a 
discrete version of AdS geometry, (b) Use of (/-invariant dis- 
entanglers and isometries, Eq. [4] guarantees that the bound- 
ary state has Q as a global symmetry, Eq. [2] (c) It also 
implies that Q acts as a local symmetry in the bulk of the ten- 
sor network (with its local action given precisely by Eq. [4|. 
Indeed, the MERA is invariant under the transformation of 
e.g. two tensors by two different group elements g and g' . 



thus obtain RG fixed points and exact MERA represen- 
tations for each element of -ff 2 (C* ,U(1)). Previously, exact 
fixed-point MERA representations were only known for 
topologically ordered phases @ . 

The above analysis is also reassuring from a numeri- 
cal perspective. Symmetric tensors are commonly used 
in MERA algorithms, because they provide several com- 
putational advantages, including faster implementations 
algorithms and access to quantum numbers [Hi]. To 
these advantages now we add the ability to identify in 
which SPT phase a given symmetric Hamiltonian H is. 
For instance, one can numerically optimize a symme- 
try protected, scale-invariant MERA (including a few 
transitional coarse-graining transformations until the RG 
fixed-point is reached) with the algorithm of Ref. [l2j ]. 
and determine the SPT phase by analyzing the resulting 
fixed-point pair (w aym ,w aym ) (23[. 

Application II: symmetryprotected holography. — The 
AdS/CFT correspondence [3j asserts the equivalence be- 
tween a conformal field theory (CFT) in d + 1 space- 
time dimensions and a theory of gravity in anti-de Sitter 
(AdS) space-time in d + 2 dimensions, where the addi- 
tional dimension in AdS corresponds to changes of scale 



4 



in the CFT. This correspondence establishes a dictio- 
nary between properties of the CFT (seen as living at 
the boundary of AdS) and properties of the gravity the- 
ory (seen to live in the bulk of AdS). In particular, the 
scaling dimensions of operators in the CFT are related 
to masses of fields in the gravity theory, and a global 
symmetry in the boundary CFT translates into a local 
symmetry in the bulk |3|. 

Being based on a RG transformation, the MERA rep- 
resents the ground state of a local lattice Hamiltonian in 
d spatial dimensions by a tensor network spanning d + 1 
dimensions, where the additional dimension also corre- 
sponds to changes in scale. [Notice that the description 
of a time-independent ground state only requires d space 
dimensions (d + 1 in holography) , instead of the d + 1 
space-time dimensions (d + 2 in holography) required to 
study time-dependent properties]. For instance, Fig.[3ta) 
shows a MERA for the ground state of a spin chain, that 
is, a lattice model in d = 1 space dimensions, and it 
spans d + 1 = 2 dimensions. This tensor network can 
efficiently describe the ground state of a quantum criti- 
cal system, which in the continuum may correspond to a 
CFT. As a matter of fact, in that case one can extract the 
conformal data of the CFT from the MERA [11]. Sugges- 
tively, the scaling dimensions of operators are obtained 
from the masses that control the exponential decay of 
correlations in the scale direction (as extracted from the 
eigenvalues of a scaling super- operator that implements 
changes of scale), which matches one of the ingredients 
of the holographic dictionary. In addition, the scaling 
of entanglement entropy in MERA was originally com- 
puted by adding contributions from different scales [HI, 
in close analogy with the holographic computation of en- 



tanglement entropy in a CFT [25[. Ref. [13| proposed 
that, indeed, the MERA should be interpreted as a lattice 
realization of the holographic principle, with the tensor 
network reproducing a discrete version of the AdS geom- 
etry. Since then, this proposal has been supported and 
extended by several authors [l4|. Here, we aim to pro- 
vide further support for this interpretation by proposing 
how to explicitly realize in the MERA another entry of 
the holographic dictionary, namely the one that relates 
global symmetries at the boundary with local symmetries 
in the bulk. 

We have argued above, see Fig. [1] that a global sym- 
metry Q of the state 1^) e y® L can be incorporated into 
the MERA by using (/-symmetric tensors. In the holo- 
graphic interpretation of the MERA, the global symme- 
try is a property of the boundary theory, see Fig.[lja),(b). 
However, using (/-symmetric tensors also implies that the 
tensor network is invariant under the action of Q on in- 
dividual tensors, where it acts simultaneously on all the 
indices of a tensor, see eq. 01 In other words, the bulk 
theory contained in the tensor network of the MERA has 
acquired a local symmetry Q, meaning that it is invari- 
ant under the action of Q on different tensors by different 



group elements, such as g and g' in Fig. [31(c). 

In most practical applications of the AdS/CFT corre- 
spondence, one exploits the fact that it is a weak-strong 
duality: when, say, the bulk theory is weakly coupled, 
and therefore can be treated within perturbation the- 
ory, then the boundary theory is strongly coupled (and 
vice- versa). It is important to realize that the MERA, 
being based only on the existence of an RG flow, of- 
fers an efficient holographic representation of most known 
ground states [26| , regardless of whether either the bulk 
or boundary theories are weakly coupled, or even whether 
the boundary theory is related to a CFT. For any such 
ground state with a global symmetry Q, we can now ob- 
tain a holographic description with an explicit local real- 
ization of the symmetry Q in the bulk. 

Discussion. — Previously, symmetric tensors had been 
used in the MERA as a means to guarantee the exact 
implementation of a global symmetry in the resulting 
many-body state. This also had other advantages, in- 
cluding the ability to target specific quantum numbers 
during a simulation, and a significant reduction in com- 
putational space and time attained by exploiting the in- 
ner structure of (/-symmetric tensors (16J. However, in 
principle there is no need to use symmetric tensors in 
order to enforce a global symmetry and, as a matter of 
fact, in some cases usin g n on-symmetric tensors results 
in a simpler description [27J- 

Here we have proposed and explored the use of sym- 
metric tensors as a means to build a symmetry protected 
renormalization group flow, and argued that this was 
necessary in order to obtain all the RG fixed-points cor- 
responding to SPT phases. We have then also shown 
that this construction produces a holographic description 
where the global symmetry at the boundary is explicitly 
realized as a local symmetry in the bulk. Entanglement 
renormalization, as a non-perturbative, real-space, lat- 
tice realization of the powerful ideas of the renormaliza- 
tion group, offers a unique testing-ground to explore var- 
ious aspects of renormalization and the holographic prin- 
ciple. In this paper, by investigating how to implement 
a global symmetry in entanglement renormalization, we 
hope to have contributed to our understanding of the role 
played by symmetries in the renormalization group flow 
and in holographic descriptions of many-body states. 
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Note: At the final stages of preparation of this 
manuscript we became aware of the recent paper by C- 
Y. Huang, X. Chen, and F.-L. Lin on "Symmetry Pro- 
tected Quantum State Renormalization" 2l|. Quantum 
state renormalization [5] and entanglement renormaliza- 
tion [H appear to be equivalent concepts. Accordingly, 
their symmetry-protected versions, as addressed in Ref. 
21] and here respectively, are also equivalent. 



5 



[1] K.G. Wilson, Rev. Mod. Phys. 47, 4, 773 (1975) M.E. 

Fischer, Rev. Mod. Phys. 70, 653 (1998). 
[2] R. Bousso, Reviews of Modern Physics 74, 825 (2002). 
[3] J. Maldace na, Adv. T heor. Math. Phys. 2, 231 (1998), 

|arXiv:hep-th/9711200| E. Witten, Adv. T heor. Math. 

Phys. 2, 253 (1998), |arXiv:hep-th/9802150| 
[4] G. Vidal, Phys. Rev. Lett. 99, 220405 (2007). G. Vidal, 

chapter in Understanding Quantum Phase Transitions, 

edited by L. D. Carr (Taylor & Francis, Boca Raton, 

2010). 

[5] X. Che n, Z.-C. Gu, X.-G . Wen, Phys. Rev. B 82, 155138 

(2010). larXlv:1103.3323l 
[6] G. Vidal, Phys. Rev. Lett. 101, 110501 (2008). 
[7] G. Evenbly and G. Vidal, Phys. Rev. Lett., 104, 187203 

(2010). 

[8] P. Corboz et al., Phys. Rev. A, 81, 010303(R) (2010). P. 

Corboz, G. Vidal, Phys. Rev. B 80, 165129 (2009) 
[9] M. Aguado and G. Vidal, Phys. Rev. Lett., 100, 070404 
(2008). R. Koenig, B.W. Reichardt, and G. Vidal, Phys. 
Rev. B, 79, 195123 (2009). 

[10] R. Koenig, E. Bilgin, Phys. Rev. B 82, 125118 (2010) R. 
N. C. Pfeifer et al. Phys. Rev. B 82, 115126 (2010) 

[11] V. Giovannetti, S. Montangero, and R. Fazio, Phys. Rev. 
Lett. 101, 180503 (20081. larXiv:0804.0520l 

[12] R. N. C. Pfeifer, G. Even bly, and G. Vidal , Phys. Rev. A 
79(4). 04030KR) (2009). larXiv:0810.0580l G. Evenbly et 
al., Phys. Rev. B 82, 161107(R) (2010). larX!v:0912. 16421 
G. Evenbly, G. Vidal, chapter in Strongly Correlated Sys- 
tems: Numerical Methods (Springer Series in Solid-State 
Sciences, 2013) 

[13] B. Swingle, P hys. Rev. D 86, 065007 (2012), 
larXiv:0905.l3T7l B. Swingle, pre-print larXiv:1209.330"4l 

[14] G. Evenbly, G. Vidal, J. Stat. Phys. 145, 891 (2011), 
larXiv: 11 06. 1082 C. Beny, New J. Phys. 15 (2013) 
023020, larXiv:1110.4872l J. Molina- Vilaplana a nd P. So- 
dano, JFEP 10, 11 (2011), larXiv:1108. 12771 H. Mat- 
sueda, M. Ishihara, and Y. Hashizume, Phys. Rev. D 87, 
066002 (2013), larXiv:1208.0206l M. No zaki, S. Ryu, and 
T. Takayanagi, JHEP 10, 193 (2012), larXiv:1208.3469l 



T. Hartman, J. Maldacena, pre-print arXiv:1303.1080 
[15] The present discussion can be generalized to states that 

transform under {V g )® L as an arbitrary linear represen- 
tation of the group Q. 
[16] S. Singh, R. N. C. Pfeife r, and G. V idal, Phys. Rev. A 

82, 050301 (2010), larXiv:0907.2"994l S. Singh, R. N. C. 

Pfeifer and G. V idal, Phys. Rev. B 83, 115125 (2011), 

larXiv: 1008.47741 S. Singh and G. Vidal, Phys. Rev. B 

86, 195114 (2012), larXiv: 1208.39191 
[17] F. Haldane, Phys. Lett. A 93, 464 (1983). I. Affleck et 

al., Phys. Rev. Lett. 59, 799 (1987). 
[18] F. Pollmann et al. Phys. Rev. B 81, 064439 (2010), 

larXiv:0910.1811l 
[19] X. Chen, Z.-C Gu, and X.-G. Wen, Phys. Rev. B 83, 

035107 (2011). larXTv:1008.3745l 
[20] N. Schuch, D. Perez-Garcia, I. Cirac, Phys. Rev. B 84, 

165139 (2011). 

[21] C.-Y. Huang, X. Chen, and F.-L. Lin, pre-print 
larXiv: 1303.41901 

[22] Strictly speaking, for each site it would be enough to 
keep a one-dimensional subspace V C (C2); ® (C2V cor- 
responding to the singlet state |V> )- 

[23] If the symmetry group Q only has a small number of 
SPT phases, one can try all possible fixed-point pairs 
(u Bym , ui Bym ) corresponding to these phases, for which we 
have described an exact expression, and restrict the nu- 
merical optimization (energy minimization) to the tran- 
sitional layers. 

[24] T his calculation w as first presented in 
|arXivxmant- ph/0610099 , which is the arxive version of 
Ref. 6]. 

[25] S. Ryu and T. T akayanagi, Phys. Rev. Lett. 96, 181602 
(2006), |arXiv:hep-th/0603001| 

[26] In some cases, such as systems in d > 2 spatial dimen- 
sions with a d — 1 dimensional Fermi surface, an effi- 
cient description may require using a generalization of 
the MERA known as branching MERA 

[27] S. Singh, G, Vidal, in preparation. 

[28] G. Evenbly, G . Vidal, pre-prints larXiv:1205.0"639l 
larXiv:1210.1895l 



